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§ : ABSTRACT 
(N 

■ Analyses of halo shapes for disc galaxies are said to give incongruous results. 

I point out that the modified dynamics (MOND) predicts for disc galaxies a 
distribution of fictitious dark matter that comprises two components: a pure 
disc and a rounder halo. The former dominates the true disc in regions of small 
^ ■ accelerations, where it controls the z-dynamics in the disc (disc flaring etc.); 

^ ■ it has a finite total mass. It also dominates the round component near the 

CN . centre where the geometry is nearly planar. The second component controls 

motions far from the plane, has a total enclosed mass that diverges linearly 
with radius, and determines the rotation curve at large radii. Its ellipticity may 
be appreciable at small radii but vanishes asymptotically. This prediction of 
Q_|. MOND differs from what one expects from galaxy-formation scenarios with dark 

Q . matter. Analyses to date, which, as they do, assume one component-usually 

^ . with a constant ellipticity, perforce give conflicting results for the best-value 

. ellipticity, depending on whether they probe the disc or the sphere, small radii 



or large ones. 

Subject headings: gravitation-galaxies: halos, kinematics and dynamics 



1. Introduction 

The modified dynamics (MOND) predicts that the global structure of the Newtonist's 
"dark halo" in disc galaxies is rather more complex than is usually assumed for such 
halos. It has a disc component and a rounder component with radius-dependent flattening, 
becoming spherical at large radii. This structure is at odds with what one naturally expects 
from existing halo-formation scenarios under the dark-matter doctrine. 

The purpose of this paper is to point out the general properties of MOND "halos" with 
the help of analytic results for simple models and situations. 
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In section 2 I explain how fictitious dark matter arises in MOND and describe its 
distribution for thin-disc galaxies. Exact expressions for the "dark-matter" distribution in 
Kuzmin discs, for which the MOND problem can be solved exactly, are given in section 
3. In section 4 I discuss generalizations of these results; in particular, the expected 
2;-distribution of the "dark matter" within the disc. Section 5 touches briefiy on implications 
for halo-shape analyses. 



2. Fictitious dark matter according to MOND 

I work with the formulation of MOND as a modification of Newtonian gravity given 



Bekenstein fc Milgrom (1984)| whereby the gravitational potential, if), of a density 



m 

distribution p* satisfies 

V-[MlWlK)W]=4vrG'p* (1) 

instead of the Poisson equation; here, Qq is the acceleration constant of MOND, p{x) ~ x 
for a; <^ 1, and p{x) ~ 1 for a; ^ 1. A Newtonist, who assumes that solves the Poisson 
equation, will deduce a mass density p = {4:nG)~^Aip, and will thus find that he requires 
"dark matter" with density 

p, = {AnG)-'A^lj-p*. (2) 

Here I investigate the distribution of this fictitious mass density for disc galaxies. 
For the sake of concreteness I consider a thin, axisymmetric, pure-disc model. The 
generalization to a galaxy model containing a bulge is straightforward and, in any case, 
does not qualitatively change the results. 

From eqs. and (Q) with p* = 0, the "dark matter" has a component outside the disc 

p^ = -{AnGy^Vin p{\Vij\/ao) ■ (3) 

It is asymptotically spherical with p^, — > {An)^^(Mao/Gy^'^R~'^ at large radii R; M is the 
total mass of the disc. This is the component that gives the asymptotically fiat rotation 
curve. There is also a thin-disc component of "dark matter" because the jump condition on 
the potential is different from the Newtonian jump for the same surface density. In MOND, 
we have across a thin layer of true surface density S* 

5[/x(g/a„)gJ = AttGE*, (4) 

where g = —Vip is the (MOND) acceleration, g^ its component perpendicular to the disc, 
and is the jump from one side of the disc to the other. A Newtonist who measure the 
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the acceleration field will deduce a surface density S = (AnG) ^^[g^]. For the ±z-symmetric 
case, which I assume all along, we have 

S = S7/i(g±/ao), (5) 

or a "dark-mater" surface density 

S,(r) = [l//i(g±K) - 1]S* = r/(r)S*(r), (6) 

where g^ is the total (MOND) acceleration just outside the disc. It follows from the 
basic assumptions of MOND that ?7 > 0, that it is very small where g 3> flo, and that 
asymptotically 77 oc r. We also see, generally, that S^, is non zero only where there is a true 
disc. To get rj^r) we first have to solve the MOND equation for g. 

Note that for a one-dimensional planar mass layer all the "dark matter" lies in the 
layer itself; outside, the MOND acceleration is constant and so gives po = 0. This means 
that near the centre of the disc, where the geometry is nearly planar, most of the dark 
matter lies in the disc and the equipotentials are highly flattened. The Newtonist's central 
surface density of the disc is 

So = u{y)^* > S*, (7) 

where y = 2'kGT,1/ ao, Sq is the true central surface density; the function z/(?/) is defined by 
z/(?/) = I~^(y)/y, where y = I{x) = (so we have viy) = l/fi[x{y)]). The argument of 

I'ly) is always of the form gu/do, where gN is some Newtonian acceleration, while x is of the 
form g/co where g is a MOND acceleration. Regardless of the exact form of fi{x) we have 
^{y) PS 1 for ?/ 1, and iy{y) ^ y~^^'^ for y -C 1. 

For the widely used form of = x(l + x^)""*^/^ we have 

K^/) = [l/2 + (r^ + l/4)^/T/^ (8) 

Equationd^ is an example where the dark- matter factor as appears e.g. in eqs.(^ 

and (^, calculated at some MOND field value x = g/ao, can be expressed in terms of 
the Newtonian-field variable g^/ao- This is possible for all one- dimensional configurations 
where an exact algebraic relation exists between the Newtonian and MOND accelerations 
gN = — Vv? and g = —Vip respectively: 

/w(g/ao)g = gN, (9) 
or inverting, the MOND acceleration is given by 

g = z/(gN/ao)gN. (10) 

I now demonstrate all the above with analytic expressions for Kuzmin discs. 
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3. Analysis of the Kuzmin disc 

Kuzmin discs have a Newtonian gravitational potential 

^ = -MG/[r' + {\z\ + hyy'\ (11) 

where M and h are the mass and scale length of the disc (we use cylindrical coordinates 
r, 2, 6). The potential above the disc {z > 0) is that of a point mass M at — h = (0, —h, 0); 
the potential below the disc is produced by the same mass at h. The surface density, Ti^^r), 
matches the jump in the 2;-gradient of the potential 



S,(r) = (2vrG)-^ 



Mh/27i{r' + = So(l + (12) 



where Eq = M/27rh'^, and u = r/h. Everywhere outside the disc the equipotential surfaces 
are concentric spheres centred at ±h. Because of this one readily gets the exact MOND 
solution ( [Brada fc Milgrom 1995|) . There the algebraic relations (P) (p!0|) hold exactly, where 



here we now have 

gj^(R) = -MG(R±h)/|R±h|^ (13) 

for the Newtonian acceleration field above (+) and below (— ) the disc. We use the 
parameter ( = MG/h^ao to measures how deep in the MOND regime the disc is. For C 1 
the whole disc is in the MOND regime, and we can write everywhere outside the disc 

^ ^ {MGaoY'nn[r' + {\z\ + hfY^\ (14) 

Otherwise this expression is only valid asymptotically (for R/h ^ C^^^)- 

The Newtonist's mass distribution comprises a disc whose surface density is 

S(r) = z/(i/)S,.(r), (15) 

where 

y{r) = g^{r, 0+)/ao = MG/a^ir' + h') = (/{l + u^). (16) 
So for Kuzmin discs 

v{r) = MC/{l+u')]-l}. (17) 

For large values of y (the Newtonian regime), as are found in the inner parts of a disc with 
C ^ 1, we have ^ 1. When -C 1, or at large radii for all (, we have 

(18) 
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and the "dark disc" dominates. (For the Kuzmin disc where S^(r) oc r ^ at large r, 
oc r^^ there, and the "dark disc" mass diverges logarithmically.) For the choice 
fi{x) = x(l +x2)-i/2 we have from eq.® r] = [1/2 + {y-^ + Xj^f^^Y'^ - 1. 

In addition we need outside the disc a rounder component of "dark-matter" with 
density 

p,(r, z) = -{AnG)-'V ■ g = -^C-'/¥/'^'(y), (19) 

where y = g^/ao = ([u"^ + (1 + ^ = and i^'{y) = du{y)/dy. We have 

Po oc asymptotically. The equi-density surfaces of this distribution coincide with the 
equipotentials of ip and ip: they are spheres centred at ±h: all quantities are functions of 
V? + {1 + vY- As expected the eccentricities are large near the centre {u,v <^ 1). 



4. Generalizations 

For a thin disc with an exponential surface density profile the algebraic relation between 
Newtonian and MOND accelerations, eqs.(|]) and (p!0|), is a very good approximation (as 
shown numerically in [Brada fc Milgrom 1995| ). Expression (|17|) for rj is thus still a good 



approximation with y{r) = g^{r,0'^)/ao] but gN is now the Newtonian acceleration for 
the exponential disc. Using the analytic expression for g.,. from [Freeman (1970)| , I find 
y{r) = ([exp{—2s) + g^(s/2)]^/^, with s = r/h, q{x) = x[Iq{x)Kq{x) — Ii{x)Ki{x)], and 
Jj, Ki modified Bessel functions. 

I have so far discussed only the z-integrated "dark-matter" surface density. We need 
to know how it is distributed in z to predict dynamics in the disc. In the limit where the 
disc is very thin the term with the radial part of the divergence in eq.(|I]) can be neglected 
in the bulk of the disc (where dgz/dz ^ dg^/dr), and we write 

p* ^ (47rG)-i^[/i(g/a„)g,]. (20) 

Or, we can turn this into an expression for the integrated surface density between ±z: 
T.*(r,z) ^ {2nG)^^ix{g/ao)gz] g = (g? + glY^'^- In our approximation g^ is independent 
of z, but gz varies strongly with z. The Newtonist's integrated surface density is 
S(r, 2;) «i {27iG)~^gz = S*(r, z)//i(g/ao). Just outside the disc gz becomes independent of z 
and we get the saturated ratio of Newtonist-to-true surface density as in eq.(|). If we are 
in a regime where g,. ^ gz for all z (this is generally valid beyond a few scale lengths) than 
g, and hence p, is approximately z independent, and the "dark matter" z distribution is 
proportional to that of the true density. Otherwise the "dark matter" is more concentrated 
near the mid-plane than the true density, since g^, and hence p, increases with z. 
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5. Discussion 



In recent papers [Sackett (1999)1 , piling fc IVlerritield (2000)1 , and llbata fc al (2000) 



note that different techniques for estimating the ellipticity of galactic dark halos give 
systematically different results. These analyses were apphed to different galaxies, so, strictly 
speaking, are not in conflict with each other. It is, however, alarming that the degree of 
flattening deduced is so strongly correlated with the type of technique used. 

But, if MOND is basically correct one cannot successfully model the dark- matter 
distribution in disc galaxies by a single- component, elliptical halo of constant ellipticity. 
Probing the galaxy far from the plane and from the centre will result in rounder halos as 
in llbata fc al (2000)| . Probing near the centre might give appreciable eccentricities. For 
example, the considerable ellipticity that ISackett fc al (1994)1 find for the halo of NGC 



4650A reflects mostly constraints from fits to the rotation curve of the smaller, main disc. 
NGC 4650A is a low- surface- density galaxy, ( 1, requiring much dark matter even at 
small radii. This, we saw, is concentrated in the disc, and so attempts to put it in a halo 
around the central disc requires a highly flattened halo. Indeed, [Vlorishima fc Saio (1995 



who analyzed the rotation curves of the main disc and the polar ring of this galaxy in a 
MOND potential derived from eq.(0) find good agreement. 

Disc-flaring analyses might also indicate a highly flattened halo. The main reason is 
that the MOND "dark disc" can greatly enhance the self gravity of the disc and make it 
dominant in determining the z structure. 

It is impossible to make quantitative comparisons with published disc-flaring analyses, 
as in IBecquaert fc Combes (1997)| and piling (1996b)l , without delving into the many 



details of the data and of the assumptions going into such analyses (e.g. the contribution of 
a stellar disc in addition to that of the studied HI layer, vertical velocity dispersion, etc.). 
As a demonstration, consider a pure disc of mass M, scale length h, and surface density 
S(r) = {M /2Tth^)f{r /h). Take r beyond a few scale lengths so that the accumulated mass 
is roughly saturated. Also assume we are in the MOND regime (^u ^ 1) and at a radius on 
the flat part of the rotation curve. Neglecting the confining effects of the inner disc and 
reckoning only with the local self gravity (i.e. assuming pure planar symmetry), the z scale 
hight according to MOND is 



where az is the z velocity dispersion. Applying instead the expression piling (1996b) 



uses 



to analyze HI flaring in NGC 4244, which neglects disc self gravity and assumes that the 
confining agent is a flattened halo, we have Zq ~ (cr^/V)[2.4g/(1.4 + q)Y^'^r (taking in his 
expression a core radius Tc <^ r), where V is the (constant) rotational velocity at r, and q is 
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the flattening. Equating the two expression we see that to get the scale hight for MOND 
with a halo one would need a q value satisfying [2.4g/(1.4 + g)]^/^ ~ {a z /V){h / rY / f {r / h) ^ 
where I made use of the MOND relation V"^ ~ MG/fir. So, for example, for an exponential 
disc, at 5 scale lengths, and with az/V = 0.1, a q value of 0.24 is gotten. 

Using data for NGC 4244 from [Oiling (1996a)| we have at the last observed point (13 
kpc) a gas surface density of S ~ 5 x 10~^gr cm~^ (corrected for helium), the MOND factor 
fi ^ V'^/rtto ~ 0.16, and the planar velocity dispersion a ~ 8km s~^. Applying the MOND 
estimate for the scale hight eq.(^), and assuming that the z dispersion equals the planar 
one, gives zq ~ 1700 pc compared with Olling's measured half width at half maximum of 
750 pc. We have to remember that cr^ might be somewhat different, that eq.(PT|) is only 
approximate, neglecting the global restoring force (see below), etc. Without the MOND 
factor ("dark disc" contribution) the expected value of zq is 6 times larger (i.e. ~ 10 kpc), 
so the true disc alone is quite negligible in containing the HI layer in Newtonian dynamics. 
At 10 kpc the MOND value of zq is ~ 300 pc (/i ~ 0.27 is larger, and S is 10 times larger, 
than at 13 kpc) compared with a similar value in [Oiling (1996a)| . So, the high degree of 



halo flattening that Oiling (1996b) flnds for this galaxy may well be largely an artifact of 



the MOND "dark disc" . 

The above was only a demonstration of the importance of the MOND effects in 
augmenting the self-gravity restoring force of the disc, even in places where the surface 
density is small. In a detailed analysis one must add the conflning force, gl, of the 
rounder "halo" (the global potential of the disc) to the local, self-gravity component Qz- 
At very large radii, or beyond some sharp cut-off of the disc, the former will dominate. 
We can estimate the relative importance of the two effects through the parameter 
f3 = Qz/ g*z ~ 27rS(r)r^/Mzo- This approximation is good inasmuch as the Newtonian radial 
acceleration can be written as MG/r^, so the MOND radial acceleration is MG/fir"^, and 
so the restoring force of a spherical potential would be zMG/^r^; also Qz ~ 27rGS(r)/yu(r). 
Applying this to the above two radii in NGC 4244, with a MOND mass of 3 x 10^ M©, 
corresponding to an asymptotic velocity of 85 km s~^, I get /? ~ 10 at 10 kpc, and f3 ~ 1.2 
at 13 kpc. This means that self gravity dominates at 10 kpc and global forces can be 
neglected, while at 13 kpc the two are of the same order, and the above MOND estimate of 
Zq at 13 kpc has to be decreased by a factor of about 2. 
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